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ABSTRACT 

We present high angular resolution observations of the Class protostar 
IRAM 04191+1522, using the Submillimeter Array (SMA). The SMA 1.3mm 
continuum images reveal within IRAM 04191+1522 two distinct sources with 
an angular separation of 7.8 + 0.2". The two continuum sources are located 
in the southeast-northwest direction, with total gas masses of ~ 0.011 M Q and 
~ 0.005 M Q , respectively. The southeastern source, associated with an infrared 
source seen in the Spitzer images, is the well-known Class protostar with a bolo- 
metric luminosity of ~O.O8L . The newly-discovered northwestern continuum 
source is not visible in the Spitzer images at wavelengths from 3.6 to 70 fxm, 
and has an extremely low bolometric luminosity (<0.03L Q ). Complementary 
IRAM N 2 H+ (1-0) data that probe the dense gas in the common envelope sug- 
gest that the two sources were formed through the rotational fragmentation of 
an elongated dense core. Furthermore, comparisons between IRAM 04191+1522 
and other protostars suggest that most cores with binary systems formed therein 
have ratios of rotational energy to gravitational energy /3 rot > 1%. This is consis- 
tent with theoretical simulations and indicates that the level of rotational energy 
in a dense core plays an important role in the fragmentation process. 

Subject headings: binaries: general - ISM: clouds - ISM: individual (IRAM 04191+1522) 
— stars: formation 



1. INTRODUCTION 

The origin of binary stars is one of the major puzzles in our understanding of star 
formation. Binary and higher-order multiple systems have been frequently observed in all 
stages of stellar evolution, but at present we do not understand well how this occurs (see, e.g., 
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reviews by Tohline 2002; Goodwin et al. 2007). The formation of binary systems begins 
in the earliest phase of the star formation process, when the young protostar, generally 
referred to as Class protostar (Andre et al. 2000), is still deeply embedded within a dense 
core of gas and dust. In order to achieve a comprehensive understanding of binary star 
formation, high angular resolution observations of the gas and optically thin dust emission 
at (sub-) millimeter wavelengths are therefore required to probe the system kinematics and 
individual core masses. However, these observations were long limited by the low angular 
resolution of single-dish millimeter telescopes, and have only become possible with the recent 
availability of large (sub-) millimeter interferometers. 

In the past decade, a handful of Class protostellar binary (protobinary) systems have 
been found using different interferometers (e.g., Looney et al. 2000; Launhardt 2004), and 
there are increasing interferometric studies of binarity in the protostellar phase (e.g., Volge- 
nau et al. 2006; Chen et al. 2007; Maury et al. 2010). Nevertheless, due to the large effort 
and observing time required, the number of known and well-studied protobinary systems is 
still small. Furthermore, high angular resolution kinematics data, which could be used to 
study the physical process leading to the formation of binary stars, are also lacking. There- 
fore, it is critical to observe, at high angular resolution, more protostellar cores in nearby 
molecular clouds to search for protobinary candidates and study in detail their kinematic 
properties, and then compare the observational results with different theoretical models. 

In this Letter, we present interferometric observations of IRAM 04191+1522 (hereafter 
IRAM 04191), using the Submillimeter Arra}0 (SMA; Ho et al. 2004). IRAM 04191 is a 
Class protostar discovered in the southern part of the Taurus molecular cloud (Andre et al. 
1999). From its low bolometric luminosity (Lboi ~ 0.15 L©) and temperature (Tb i~18K), 
Andre et al. (1999) estimated an age of ~ 1-3 x 10 4 yr, suggesting IRAM 04191 as one of the 
youngest protostars in Taurus. IRAM 04191 drives a large scale (~0.1pc), bipolar molecular 
outflow, which extends in the northeast-southwest direction (Andre et al. 1999; Lee et 
al. 2005). Various molecular line observations have shown subsonic infall, fast rotation, 
depletion, and deuteration in IRAM 04191 (see, e.g., Belloche et al. 2002; Takakuwa et al. 
2003; Belloche & Andre 2004; Lee et al. 2005). This source is also detected in the Spitzer 
Space Telescope (Spitzer) images at wavelengths between 3.6 and 70 /im, and is classified as 
a very low luminosity object (VeLLO), due to its low internal luminosit}@ (L int ~ 0.08 L Q ; 



1 The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and 
the Academia Sinica. 

2 The internal luminosity is the luminosity of the central source, which excludes the luminosity arising 
from external heating. 
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Dunham et al. 2006, 2008). Based on the SMA observations, we report in this Letter that 
IRAM 04191, regarded as an excellent example of isolated single star formation in the past 
decade, is actually a binary system, which may in turn represent one of the best examples 
of a protobinary system formed though rotational fragmentation. 

2. OBSERVATIONS AND DATA REDUCTION 

The SMA 230 GHz observations of IRAM 04191 were carried out in 2007 November (~ 5 
hours integration time) with the compact configuration (PI: J. Karr), and in 2011 March (~ 1 
hour integration time) with the subcompact configuration (PI: X. Chen). Seven antennas 
were used in both configurations, providing baselines ranging from 9-53 kA (compact) and 6- 
51 kA (subcompact), respectively. In the 2007 observations, the digital correlator was set up 
to cover the frequencies from 219.5 to 221.4 GHz and from 229.5 to 231.4 GHz in the lower 
and upper sidebands (LSB and USB), respectively. The 1.3mm continuum emission was 
recorded with a total bandwidth of ~ 3.12 GHz, combining the line-free portions of the two 
sidebands. In the 2011 observations, the SMA bandwidth was upgraded to 4 GHz, and the 
correlator was set up to cover the frequencies from 216.8 to 220.8 GHz (LSB) and from 228.8 
to 232.8 GHz (USB), respectively, resulting in a total continuum bandwidth of ~ 7.46 GHz. 
The SMA primary beam is about 55" at 230 GHz. 

The visibility data were calibrated with the MIR package (Qi 2005). In the 2007 ob- 
servations, the quasar 3c454.3 was used for bandpass calibration, and quasars 3cl20 and 
0530+135 for gain calibration. In the 2011 observations, the quasar 0854+201 was used for 
bandpass calibration, and quasars 3c84 and 0359+509 for gain calibration. Uranus was used 
for absolute flux calibration in the 2007 observations, while 0854+201 was used in the 2011 
observations. We estimate a flux accuracy of ~20% for the data, by comparing the final 
quasar fluxes with the SMA calibration database. The calibrated visibility data were further 
imaged using the Miriad toolbox (Sault et al. 1995). 

3. RESULTS 

Figure 1 shows the SMA 1.3 mm dust continuum images (top panel) and the Spitzer 
infrared images (bottom panel) of IRAM 04191. The SMA 1.3 mm dust continuum observa- 
tions reveal within IRAM 04191 two distinct sources, which are separated by 7.8 + 0.2" in 
the northwest-southeast direction (see Figures la-c). These two sources were detected with 
the SMA compact and subcompact configurations independently, thus yielding a solid de- 
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tection. The SMA compact configuration observations detected relatively compact emission 
clearly arising from two distinct sources (see Figure la), while the subcompact configuration 
observations detected relatively extended emission with two local maxima (see Figure lb), 
likely tracing an inner common envelope. As shown in the Spitzer images (see Figures ld- 
f), the southeastern continuum source is associated with an infrared source, and is therefore 
referred to as IRAM 04191 IRS in this work (see also Dunham et al. 2006). In contrast, the 
northwestern continuum source has no compact infrared emission at the Spitzer bands from 
3.6 to 70 jum, and is referred to as IRAM 04191 MM£§|. The comparisons with previous high 
angular resolution observations (e.g., Belloche et al. 2002; Lee et al. 2005) suggest that 
source IRS is the well-known Class protostar, which drives a large scale bipolar outflow in 
the northeast-southwest direction. 

The continuum fluxes of the two sources were derived from Gaussian fits to the restored 
images (using the combined compact and subcompact data) with the Miriad command imfit 
(see Table 1). The fluxes of sources IRAM 04191 IRS and MMS are ~ 9.3 mJy and ~ 4.7mJy, 
respectively. For comparison, the continuum flux detected in the IRAM-30m 1.3 mm map is 
~650mJy (aperture 60"; Andre et al. 1999), indicating that more than 95% of the flux from 
the large-scale envelope around the two sources was resolved out by the SMA observations. 
Furthermore, it must be noted that in the IRAM Plateau de Bure Interferometer (PdBI) 
1.3 mm continuum observations, only a single object was detected in IRAM 04191, which 
is spatially coincident with source IRAM 04191 IRS (see Maury et al. 2010). It is likely 
that the PdBI observations missed the relatively faint source MMS (with about half of the 
IRS flux), due to the low flux ratio sensitivity in the IRAM 04191 observations (only about 
0.6-0.9; Maury et al. 2010). 

Assuming that the 1.3 mm dust continuum emission is optically thin, the total gas mass 
(M gas ) of the two sources is derived from their fluxes using the same method as described 
in Launhardt & Henning (1997). In the calculations, we adopt a dust opacity of = 
0.5 cm 2 g -1 (Ossenkopf & Henning 1994), a typical value for dense and cold molecular cloud 
cores. A dust temperature of ~20K was adopted for IRAM 04191 MMS and IRS, which is 
derived from fitting the SEDs (see below). The relative uncertainties of the derived masses 
due to the calibration errors of the fluxes are ±20%. The total gas masses of IRAM 04191 IRS 
and MMS derived from the SMA maps are ~ 0.011 M and ~ 0.005 M , respectively. 



3 The Spitzer spatial resolution is ~ 18" at 70 /xm, which fails to resolve sources IRS and MMS. However, 
due to the facts that (1) the 70 /xm emission peaks at source IRS, without any shift toward source MMS, 
and (2) the emission is roughly circularly symmetric (see Figure If), we consider that all the 70 fim emission 
comes only from source IRS. 
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4. DISCUSSION 

4.1. The Spectral Energy Distributions 

In order to understand the properties of sources IRAM 04191 IRS and MMS, it is es- 
sential to construct their spectral energy distributions (SEDs), and then to derive their 
luminosities and temperatures. However, previous single-dish (sub-) millimeter continuum 
observations lack the resolution to resolve the two sources. Therefore, we simply assume 
a flux ratio of 2:1 at these wavelengths (from 160 /xm to 1.3 mm), which is inferred from 
the SMA flux ratio of the two sources. Furthermore, due to the fact that no compact in- 
frared emission was detected from MMS in the Spitzer 3.6 to 70 ^m bands but infrared 
emission was detected at the position of IRS, we assume that all the infrared flux at wave- 
lengths <90fim comes from source IRS. Given these assumptions, we show in Figure 2 the 
SEDs of sources IRS and MMS (see Dunham et al. 2008 for individual flux values). To 
derive luminosities and temperatures, we first interpolated and then integrated the SEDs, 
assuming spherical symmetry. Interpolation between the flux densities was done by a \ 2 
single-temperature grey-body fit to all points at A > 90 ^m, using the method described in 
Chen et al. (2008). A simple logarithmic interpolation was performed between all points at 
A <90 fim (including 3 a upper limits for source MMS). From the SED fitting, we estimate 
^boi ~ 0.08 L Q and T bol ~ 25 K for source IRS, and L ho \ < 0.03 L Q and Tboi < 20 K for source 
MMS. The dust temperature estimated for the two sources is 20±3K. Further observations, 
such as high angular resolution and high sensitivity continuum observations at wavelengths 
from far-infrared to (sub-) millimeter, are needed to constrain the SEDs of IRAM 04191 IRS 
and MMS, in order to derive more precisely their luminosities and temperatures. 



4.2. A Case of Rotational Fragmentation 

It has been suggested that the fragmentation of rotationally flattened cloud cores with 
initially flat density profiles, immediately after a phase of free-fall collapse, is one of the most 
efficient mechanisms for binary star formation (see reviews by Tohline 2002 and Goodwin 
et al. 2007). Numerical simulations find that rotating cloud cores fragment if a t herm/3 r ot 
< 0.12-0.15 (e.g., Tohline 1981), where atherm = -E'therm/-E , grav is the initial thermal virial 
ratio (.Etherm is the thermal kinetic and E gra , v is the gravitational potential energy) and /3 rot 
= E TOt / E gTay is the initial rotational virial ratio (E TOt is the rotational kinetic energy). In 
particular, the simulations find that the level of rotation energy, i.e., /3 rot , plays an important 
role in this fragmentation process. For example, Boss (1999) found that rotating cloud cores 
fragment only when /3 rot > 0.01 (see also Machida et al. 2005). 
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Figure 3 shows the velocity field of IRAM 04191, using the combined IRAM-30m and 
PdBI N 2 H + (1-0) data from Belloche & Andre (2004) and derived with the same method 
used by Chen et al. (2007). The mean velocity map shows a clear continuous velocity 
gradient across the dense core, increasing from northwest to southeast, i.e., roughly along 
the line connecting sources MMS and IRS (see Figure 3a). As discussed in Chen et al. (2007), 
systematic velocity gradients found in protostellar cores/envelopes are usually dominated by 
either rotation or outflow. As seen in Figure 3a, the large angle (> 70°) between the gradient 
and the axis of the IRS outflow suggests that the observed velocity gradient in N 2 H + is due 
to core rotation rather than outflow. Figure 3b shows the spatial distribution of the N 2 H + 
line widths in the map. The line widths are roughly constant within the core, and relatively 
large line widths are mainly seen in the central region of the core (roughly along the outflow 
axis), which might be caused by the outflow (see Chen et al. 2007 for more discussions). 

From the velocity field derived from the N 2 H + (1-0) observations, we can estimate the 
ratios of rotational and thermal energy to gravitational potential energy, using equations 
obtained from Chen et al. (2007; 2008): 

** = 0263 x # 2 x x &• « 

tttherm = 0.01 X (-) X )\ (2) 

K Av ohs 

where R is the core radius (the FWHM radius ~20" or 2800 AU), g is the velocity gradient 
(~ 17kms~ 1 pc -1 ), T is the kinetic gas temperature (6-10 K; see Belloche et al. 2002), and 
Af bs is the observed mean line width (0.50-0.55 kms -1 ; see Figure 3b). The estimated 
/3 TOt and a t hcrm values for IRAM 04191 are 0.06 ±0.02 and 0.3 ±0.1, respectively. Correcting 
for inclination (' = 50°; see Belloche et al. 2002), the estimated /3 ro t is then about 0.17, 
consistent with the value of ~ 0.1-0.2 derived by Belloche et al. (2002). It is clear from 
these values that IRAM 04191 meets the classical fragmentation criterion of a t hcrm/3rot < 
0.12, and the estimated (3 rot value is far beyond the critical value suggested by numerical 
simulations. We thus suggest that the two sources in IRAM 04191 were formed through 
the rotational fragmentation of an elongated dense core. Furthermore, because of its clear 
rotational picture, close distance, and youth, as well as its other well-characterized properties 
derived from previous observations, IRAM 04191 represents so far one of the best examples 
to observationally study the rotational fragmentation of dense cores. 

Interestingly, Tohline (1981) pointed out that the semi-major axis of a binary formed 
through rotational fragmentation is 

4 

a bin — ( — )/5rot-Rcorc- (3) 
7T 
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For IRAM 04191, its core radius derived from the IRAM-30m 1.3 mm observations is about 
4200 AU or 30" (see Andre et al. 1999). Adopting f3 TOt = 0.17, the predicted binary separation 
is ~6.5", close to the value measured in the SMA 1.3mm continuum images (~7.8"). 

4.3. The Role of Rotation in Binary Star Formation 

At this stage, it is instructive to compare IRAM 04191 to other protostars with high 
angular resolution kinematics data, and then to observationally investigate the role of ro- 
tation in binary star formation. Based on our previous observations (Chen et al. 2007; 
2008; 2009), we have collected high angular resolution N 2 H + (1-0) data for twelve Class 
protostars. Among those protostars, two protobinary systems, CB 230 (Chen et al. 2007) 
and SVS 13B (Chen et al. 2009), show a similar rotational picture as seen in IRAM 04191, 
and may represent another two cases of rotational fragmentation. In addition, VLA NH3 
observations suggest core/envelope rotation in the protostars HH212-MMS (Wiseman et al. 
2001) and HH211-mm (Tanner & Arce 2011), both of which were also resolved into binary 
systems in high angular resolution observations (Codella et al. 2007; Lee et al. 2009). More 
recently, Tobin et al. (2011) present high angular resolution N 2 H + and NH 3 data for another 
eleven protostars. With these kinematic data, we estimate the /3 rot values for the individual 
sources (using the same method given above), and show in Figure 4 the distribution of /3 TOt 
versus core FWHM radius for those sources (more analyses based on these kinematics data 
will be presented in a forthcoming paper). 

As seen in Figure 4, we find that most cores with binary systems formed therein have 
/3 rot > 1%. In contrast, those cores with /3 rot < 1%, such as the well-studied protostar 
LI 157, remain as single protostars in the observations. This is consistent with the theoretical 
simulations (see above) and indicates that the level of rotational energy in a dense core plays 
an important role in the fragmentation process. Nevertheless, a few protostars with /3 TOt > 
1% also remain as single (see Figure 4). We speculate that these protostars are unresolved 
binary/multiple systems due to the limitation of angular resolution or sensitivity. Indeed, 
for example, our SMA HCO + and CO observations suggest that one such protostar, L1521F, 
is likely a multiple system with a quadrupolar outflow (X. Chen et al. in preparation). 
However, it must be noted that (1) the size of the sample in Figure 4 is still small and no 
statistically significant conclusion can be drawn yet; and (2) more high angular resolution 
observations are needed to study a larger sample thoroughly and to investigate if other core 
properties, e.g., turbulence, density structure, and magnetic field, also have effects on the 
fragmentation process. 
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5. SUMMARY 

We report the discovery of a binary system with an angular separation of 7.8 ±0.2" in 
the well-known Class protostar IRAM 04191, using the SMA observations. Complementary 
IRAM N 2 H + (1-0) data suggest that the binary system is formed through the rotational 
fragmentation of an elongated dense core. Furthermore, comparisons between IRAM 04191 
and other protostellar cores suggest that most cores with a binary system formed therein 
have ratios of rotational energy to gravitational energy f3 TOt > 1%, which is consistent with 
theoretical simulations and indicates that the level of rotational energy in a dense core plays 
an important role in the fragmentation process. 

We thank the anonymous referee for many insightful comments and suggestions. We 
thank A. Belloche and P. Andre for providing their N 2 H + (1-0) data of IRAM 04191 + 1522. 
This material is based on work supported by NSF grant AST-0845619 to HGA. 
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Table 1. SMA 1.3 mm dust continuum results of IRAM 04191 



Source 


R.A. & Dec. (J2000) 


Flux a 


FWHM Sizes" 




Name 


[h : m : s, ° : ' : "] 


[mjy] 


Maj. X Min. 


[Mq] 


IRAM 04191 IRS 
IRAM 04191 MMS 


04:21:56.90, 15:29:46.4 
04:21:56.37, 15:29:48.9 


9.3±2.0 
4.7±1.0 


3.6" x 3.5" 
2.5" x 1.6" 


0.011±0.002 
0.005±0.001 



a Flux and FWHM sizes derived from the Gaussian fitting. 
b Total gas mass; See text for dust temperature and opacity used. 
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Fig. 1. — (a) SMA 1.3 mm dust continuum image of IRAM 04191, taken with the SMA 
compact configuration. The contours start at ±3 a, and then increase from +3 a in steps of 
~1<7 (1 (j~0.70 mJybeam -1 ). The synthesized SMA beam is shown as a grey oval in the 
bottom right corner, (b) Same as Figure a, but taken with the SMA subcompact configu- 
ration (la^l.emJybeam" 1 ). (c) The SMA 1.3 mm continuum image, using the combined 
compact and subcompact data (1 a ~ 0.77 mJy beam -1 ), (d) The Spitzer IRAC 8.0 /xm im- 
age of IRAM 04191. The two crosses show the positions of the two SMA continuum sources 
IRS (white) and MMS (black), (e) Same as Figure d, but the MIPS 24 fim image, (f) Same 
as Figure d, but the MIPS 70 tun image (the contours start at 16MJysr _1 , and increase in 
steps of lOMJysr -1 ). The Spitzer data of IRAM 04191 are adopted from Dunham et al. 
(2006). 
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Fig. 2.— Spectral energy distributions of IRAM 04191 IRS and MMS. Solid lines show the 
best-fit for all points at A > 90 /xm using a grey-body model S u = B u (T d )(l — e~ T ")il, where 
B u (Td) is the Planck function at frequency v and dust temperature Td, r v is the dust optical 
depth as a function of frequency r oc u 18 , and Q is the solid angle of the source. Dashed lines 
show a simple logarithmic interpolation performed between all points at A < 90 fim for the 
two sources (the datapoints of MMS at wavelengths < 90 yum are 3 a upper limits derived 
from the infrared images). 
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Fig. 3.— N 2 H+(l-0) velocity fields of IRAM 04191+1522, derived from the N 2 H+ data 
combining the IRAM-30m and PdBI observations (data from Belloche & Andre 2004). The 
top map shows the mean velocity field (color shades), while the bottom map shows the 
spatial distribution of line widths (color shades). In both maps, black contours show the 
integrated intensity, from 0.16 to 2.1 in steps of 0.16 Jybeam -1 kms -1 (~6 times the rms 
noise); white crosses show the positions of the SMA continuum sources IRS and MMS; red 
and blue arrows show the directions of the red- and blue-shifted CO outflow lobes. The 
synthesized cleaned beam in the N2H + observations is shown in the bottom right corner. 
In the mean velocity map, the black arrow in the bottom right shows the direction of the 
velocity gradient (P.A. ~133°). 
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Fig. 4. — Ratio of rotational energy to gravitational energy /3 rot vs. core FWHM radius. 
Data for the sources shown are adopted from Chen et al. (2007; 2008; 2009) and Tobin et al. 
(2011), except IRAM 04191 (this work), NGC 1333 IRAS 4 A (Belloche et al. 2006), HH212- 
MMS (Wiseman et al. 2001), HH211-mms (Tanner & Arce 2011), and L1448C/L1448 IRS3 
(X. Chen et al. in preparation). Note that except for IRAM 04191, the values of f3 rot in this 
plot are lower limits, as they are not corrected for inclination. 



